Effect of thermal fluctuations on spin degrees of freedom in spinor Bose-Einstein 

condensates 
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We consider the effect of thermal fluctuations on rotating spinor F — 1 condensates in axially- 
symmetric vortex phases, when all the three hyperfine states are populated. We show that the 
relative phase among different components of the order parameter can fluctuate strongly due to the 
weakness of the interaction in the spin channel. These fluctuations can be significant even at low 
temperatures. Fluctuations of relative phase lead to significant fluctuations of the local transverse 
magnetization of the condensate. We demonstrate that these fluctuations are much more pronounced 
for the antiferromagnetic state than for the ferromagnetic one. 

PACS numbers: 03.75.Kk, 03.75.Lm 
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I. INTRODUCTION 



Properties of rotating spinor Bose-Einstein conden- 
sates attract a lot of attention now. First examples of 
these systems with hyperfine spin F = 1 were found 
in optically trapped 23 Na 0]. Vortex phase diagram of 
spinor condensates is very rich, since the order parameter 
has three components in F = 1 case and five components 
in F = 2 case. Topological excitations in spinor con- 
densates were studied theoretically in a large number of 
articles see, e.g., Refs HUH HQ. 

At the same time, an interest is now growing to tem- 
erature effects in atomic condensates. Refs. 0, H, 
Hfl Hll | study theoretically the Berezinskii-Kosterlitz- 
Thouless (BKT) transition associated with the prolifer- 
ation of thermally-excited vortex- antivortex pairs. For 
instance, in Ref. [a] it was shown that in quasi two- 
dimensional condensates BKT transition can occur at 
rather low temperatures, T ~ 0.5T C , at number of parti- 
cles in the system N ~ 10 4 . Recently, some signatures of 
possible BKT phase were also found close to the critical 
temperature T c in experimental work |12| . where con- 
densates in optical lattice have been studied. Finally, 
experimental evidence for the BKT transition in trapped 
condensates was reported in Ref. |l|$. Refs. pHIT^ deal 
with the thermal fluctuations of positions of vortices in 
rotated scalar condensates. Note that, according to the 
Mcrmin-Wagner-Hohenberg theorem, Bose-Einstein con- 
densation is not possible in 2D homogeneous systems. 
However, application of the trapping potential leads to 
the macroscopic occupation of the ground state of Bose 
gas. 

The aim of the present paper is to study the effect 
of thermal fluctuations in rotated quasi two-dimensional 
spinor condensates. These systems have a specific degree 
of freedom, associated with the relative angle among dif- 
ferent components of the order parameter corresponding 
to different hyperfine state. In other words, this angle de- 
termines coherence among components of the order pa- 
rameter. Also it influences a transverse magnetization 
of the condensate. In this paper, we focus on thermal 



fluctuations of this angle. Note that experimentally, at 
present time, it is possible to study the condensate phase 
|l6l fill IT^ . see also Ref In addition, recently, a 

new and nondestructive method for measuring the local 
magnetization of the condensate was proposed and suc- 
cessfully applied in Ref. ppf . 

We show that the relative angle among hyperfine com- 
ponents of the order parameter in 2D case can expe- 
rience strong thermal fluctuation even at low tempera- 
tures. The reason is the weakness of the spin energy of 
the system as compared to interactions in density chan- 
nel. Also fluctuations of this angle lead to significant rel- 
ative fluctuations of the local transverse magnetization 
of the condensate, which are much larger in the antifer- 
romagnetic case than in the ferromagnetic one. 

This paper is organized as follows. In Section II, we 
give a basic formulation of the problem. In Section II, 
we discuss our main results for the fluctuations of angle 
and spin textures. We conclude in Section III. 



II. BASIC FORMULATION 

We consider harmonically-trapped quasi 2D Bose- 
Einstein condensate with spin F = 1. The trapping po- 
tential is given by 



U(r) 
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where uj± is a trapping frequency, m is the mass of the 
atom, and r is the radial coordinate. The system is ro- 
tated with the angular velocity fl, well below the critical 
rotation speed lu± , and the number of atoms in the cloud 
is N. In this paper, we restrict ourselves on the range of 
temperatures much smaller than T c . Therefore, we can 
neglect a noncondensate contribution to the free energy 
of the cloud. The total energy of the system in this ap- 
proximation coincides with the energy of the condensate. 
For the number of condensed particles, we use the ideal 



2 



gas result: 



N(T) = N 



1- — 



At the same time, 



the same, and the themperature to fix dimensionless t. 
In real atomic condensates, a is approximately several 
nanometers. The most realistic value of N for this g n 
>■ > is close to 10 3 , and to illustrate the effect of N we will 
consider the following range: 10 2 < N < 10 4 . 

The total magnetization of the condensate is fixed: 



Hu)± 
kBT c 



cm 
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where £(2) is a Riemann zeta function, yC(2) « 1.28. 
Eqs. (2) and (3) remain accurate even for the case of 
interacting particles |2l|. We also introduce a dimen- 
sionless temperature t = T/T c . Since we are considering 
low temperatures, T ~ 0.1T C , temperature dependence of 
condensed particles number can be neglected, N(T) ~ N. 

The order parameter in the F = 1 condensate has three 
components tyj (j = —1, 0, 1). The free energy of the 
system can be written as [23, |23| 



F = hu ± N / dS 



2TTg s ^* l {F a ) 0k {F a )i m ^ k ^ m 

•*n-«*(v x r)%] , 



(4) 



where the integration is performed over the system area, 
repeated indices are summed, F a (a = x, y, z) is the an- 
gular momentum operator, which can be expressed in a 
matrix form through the usual Pauli matrices, h is the 
one-body Hamiltonian, given by 



h 
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Constants g n and g s characterize interactions in den- 
sity and spin channels and are given by 



9s 



(a + 2a 2 )n 2 



(02 - a )n z 



(6) 



(7) 



where ao and are scattering lengths for atoms with to- 
tal spin and 2, and n z is the concentration of atoms in 
longitudinal direction. In real spinor condensates, \g„\ -C 
|<7n|) since ao ~ a%. Typically, \g s /gn\ ~ 0.001 — 0.01, and 
this ratio can be tuned. In this paper, we study the case 
of relatively dilute condensate and take g n = 10. We will 
consider different values of N but at fixed value of inter- 
action parameter g n . This is possible, since, in the case 
of a single layer cloud, we can always tune the trapping 
frequency in the longitudinal direction keeping g n con- 
stant. To ensure the regime of quasi-two-dimensionality, 
we can also tune lo± . In this case, we have to change the 
rotation speed to keep the dimensionless rotation speed 



M 



dS\9j\ 2 j. 



(8) 



Magnetization M is normalized in terms of N and max- 
imum of \M\ is equal to 1. One has also to take into 
account the normalization condition for the order param- 
eter: 



(9) 



The spatial profiles of all the components of the order 
parameter in the equilibrium can be found from the con- 
dition of minimum of energy (4). It is also convenient 
to introduce the longitudinal l z and transverse lt r local 
magnetizations of the condensate: 



l*i| 



(10) 
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= 2|* | 2 |* 1 | 2 + 2|^o| 2 |*_ 1 | 2 



+4 x + C.C.) . (11) 

The spin energy in this case can be represented as 



F sp in = 27r g s hjj±N J dS [l\ + I 



(12) 



In this paper, we restrict ourselves only to the case of 
axially-symmetric phases, when moduli of all the com- 
ponents of the order parameter are independent on the 
azimuthal angle and depend only on the radial coordi- 
nate r. Note that equilibrium vortex phases in this sit- 
uation were studied in Refs. 0, for the spin F = 1 
condensate and in Ref. 0] for F — 2 system. For the 
axially-symmetric phases, each component of the order 
parameter can be represented as 



*j ( r ) <P) = fj( r ) exp(-iLj<p - i5j), 



(13) 



where ip is a polar angle, Lj is a winding number, and Sj 
is a relative phase. We will denote such phases as (L-i, 
Lq, Li). As it was shown in Ref. 0, an axial symmetry 
of the solution implies that winding numbers satisfy the 
relation: L\ = ILq. In this case, according to Eqs. 

(11) and (12), the spin energy depends on relative angle 
X = 2S - Si - 6-i. 
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It is important to note that only a spin contribution to 
the total energy (4) depends on phases Sj via the spin- 
mixing term. For the stationary state, which is a local 
minimum of Gross-Pitaevskii functional (4), the value of 
X is determined by the sign of interaction constant in a 
spin channel g s . For positive g s (antiferromagnetic case), 
a minimum of F sp i n is attained at X = whereas for 
negative g s (ferromagnetic case) x = 0. 



III. RESULTS AND DISCUSSION 

According to the results of Ref. @, for the antifer- 
romagnetic state (g s > 0), phases (-1, 0, 1) and (1, 1, 
1) are energetically favorable in the region of small and 
moderate values of magnetization M. Phase (-1, 0, 1) 
is realized at low rotation frequencies f2 and (1, 1, 1) at 
higher f2. In Ref. 0, it was shown that (0, 1, 2) state 
is favorable in the ferromagnetic case (g s < 0) in a re- 
gion of moderate values of fl and M. In these phases, 
all three hypefme states are populated. Fluctuations of 
X have a sense only in this case, since the ^-dependent 
part of the energy is equal to zero identically, if one of 
the components of the order parameter is zero. In this 
paper, we will concentrate on these three vortex states, 
since they are appropriate candidates for the illustration 
of the effect of thermal fluctuations. Note that in ho- 
mogeneous spin-1 condensate atoms populate only two 
or one hyperfme state(s); they can populate three states 
only if the system is trapped and experiences rotation, 
which generates vortices. 

An important feature of real atomic spinor Bose- 
Einstein condensates is a weakness of the spin interac- 
tions comparing to the interaction in density channel 
( | .g s | <C \g n \)- At the same time, the coherence among 
the different components of the order parameter (angle 
X) is fully determined by the spin interaction. Angle x 
also influences the transverse magnetization of the con- 
densate, as seen from Eq. (11). Note that a longitudinal 
component of magnetization is independent on x- 

Smallncss of g s comparing to g n leads to the fact that 
thermal fluctuations of relative angle x become signif- 
icant at much lower temperatures than fluctuations of 
the density of particles. Therefore, at relatively low tem- 
peratures, one can assume that the moduli of all the 
components of the order parameter remain fixed (that 
can be also checked numerically), whereas x is fluctu- 
ating. For the case of small fluctuations of x, one can 
use a harmonic approximation and represent the devia- 
tion of the energy of the system from the equilibrium, 
SF = F(xo + Sx) — F(xo), as a quadratic function in 
terms of the deviation of angle x from the equilibrium 
<5x = X - Xo: 



SF = 2wg s flui±N I (cos(xo + Sx)-cos(xo)) 
« 7r\g s \hio ± NI(S X f, 



where / = J dS (fif-if§ ) ■ Under these assumptions, 
the average square of the deviation of x from the equi- 
librium is given by 



((Sx) 2 ) T = 



Jd(S x ) (*x) 2 exp(-^) 
Jd(S x )eM-^r) 



(15) 



Integrals in Eq. (15) can be calculated analytically. After 
taking into account Eq. (3), we get: 



((SX) 2 ) T = 



t 



1.28VW4 & I 



(16) 



We also introduce a quantity Ax = \J ((Sx) 2 )t> which 
can be considered as an average deviation of angle X from 
the equilibrium. We see that Ax depends on dimension- 
less temperature t — T/T c , the number of particles N 
and integral /. For a given vortex phase, / is also a 
function of total magnetization M. It is important to 
emphasize that the scaling relation (16) has a sense only 
if g n is independent of N, as discussed above. In order 
to calculate /, we use a variational method, which was 
previously applied by us in Ref. Q to evaluate energies 
of various axially-symmetric vortex phases in spin F = 2 
condensate. In this approach, each component of the or- 
der parameter is modeled by a trial function and values 
of variational parameters are found from the condition of 
minimum of total energy. 

In Fig. 1 we plot calculated dependence of Ax (mea- 
sured in degrees) as a function of the number of particles 
in the system for different vortex phases at t = 0.1 and 
g n = 10. This value of g n is close to typical experimental 
ones (ao ss 5 nm, n z w 2 nm -1 ), see also calculations 
of Ref. 0,0. We assume that, for (-1, 0, 1) and (1, 
1, 1) states, g s — 0.01.g„ and M — 0.1, whereas for (0, 
1, 2), g s = -0.01p„ and M = 0.5. Note that Ax for 
the particular phase is independent on f2, since / has the 
same property. We see that even for quite low tempera- 
tures, Ax can be rather large and the coherence among 
different components of the order parameter is practi- 
cally destroyed. For smaller value of \g s \, fluctuations of 
X are, of course, even stronger. To illustrate the effect 
of temperature, in the inset to Fig. 1, we show the de- 
pendence of Ax on T for (0, 1, 2)-phase (curve 1), (1, 
1, 1) phase (curve 2), and (-1, 0, 1) phase (curve 3) at 
fixed number of atoms N — 1000, g s — —0.05g n for the 
first curve and g s = 0.05g n for two others. Note that Ax 
is almost independent on total magnetization M of the 
condensate. 

As we already pointed out, fluctuations of x lead to 
that of Itr- In a harmonic approximation, one can ex- 
press the average deviation of \lt r \ from the equilibrium 
(5 1 lt r | )t through the deviation of x : 



(S\hr\h 



(-ir.(fef), 



hf-i 



(14) 



(A + (-i)" +1 /- 



(17) 
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where u = for the antiferromagnetic case and u = 1 
for the ferromagnetic one. If in the antiferromagnetic 
state the total magnetization is not large, M < 0.5, one 
can expect that — /-i) 2 <C fif-i, and, therefore, even 
small fluctuations of \ lead to strong relative fluctuations 
of |Ztr|. At the same time, for ferromagnetic case, u = 1 in 
this equation, and relative fluctuations of \lt r \ are much 
smaller. 

We have calculated (S \hr\)T f° r different vortex phases 
and our calculations revealed that (S \ltr\)x / \hr\ is al- 
most independent on radial coordinate r for vortex 
phases (-1, 0, 1) and (1, 1, 1). This is due to the 
fact that = \Li\ for these states, therefore, fi(r) 

is nearly proportional to /_i(r), and, according to Eq. 
(15), (S \ltr\) T I \hr \ should only slightly depend on r. In 
Fig. 2 we present (5 \lt r \) T / \hr\ as a function of to- 
tal magnetization of the condensate for (-1, 0, 1) state 
at t = 0.1, g s — 0.01g n (antiferromagnetic case), and 
N = 1000. We see that relative fluctuations of trans- 
verse magnetization can be significant even at low tem- 
perature. Value of (8 \lt r \)x I \hr\ decreases with increase 
of M. This result is natural, since condensate becomes 
more polarized with growing M. An absolute value of 
((5|Zt r |) T also remains sizable. Although value of frac- 
tional quantity (5 |Ztr|)y / \hr\ is growing with decreasing 
of M, the value of \lt r \ itself becomes smaller. Therefore, 
we found that the most appropriate values of M to ob- 
serve fluctuations of transverse magnetization is around 
AI = 0.2, where both (S IhrDr I \^tr\ and |Zt r | are high: 
(6 \l tr \) T / \kr\ ^ 0.1, whereas l tr is comparable to the 
longitudinal magnetization l z in the fully polarized state 
at M = 1, where it should be easily detectable exper- 
imentally. Value of (5 \ltr\) t I \hr\ depends also on the 
vortex phase; we found that in (1, 1, 1) state it is even 
much larger than in (-1 ,0 ,1) state. 

Also we have calculated (8 \ltr\)-r for the ferromagnetic 
(0, 1, 2) phase. As can be expected, in this case, relative 
fluctuations of \lt r \ are much weaker. Physically, this is 
because \lt r \ is proportional to the ferromagnetic order 
parameter which is responsible for the ferromagnetic 
ordering. Therefore, one can expect that in the ferro- 
magnetic phase this order parameter is more robust with 
respect to thermal fluctuations, than in the antiferromag- 
netic one. In addition, average deviation of \l tr \ from the 
equilibrium is negative and its modulus is growing with 
increase of M, in contrast to the antiferromagnetic sys- 
tem. 

Thermal fluctuations should also be important in the 
case of F = 2 condensate, where there are two interaction 
constants in spin channel and two characteristic angles. 
Therefore, one can expect more complicated behavior, as 
compared to F = 1 condensate. For instance, in homoge- 
neous F = 2 system, a cyclic state can have a lowest en- 
ergy; in this case atoms populate three hyperfine states, 
and the spin energy depends on the coherence among 
them. Fluctuation problem for this system was analyzed 
in Ref. |24|. A new method to create such entangled 
states in spin-1 condensate was recently applied exper- 



imentally in Ref. |19j, where a microwave energy was 
injected to the system. As a result, particles redistribute 
from spin —1 state to spin and 1 states, and all three 
magnetic sublevels become populated. The spin-mixing 
dynamics in F = 1 condensate was studied theoretically 
in Ref. 

Note that in Eq. (14) we have assumed that fluctuat- 
ing x is spatially independent that is not true in general 
case. However, spatial gradients of x gi ye some addi- 
tional contribution to the kinetic energy of the system, 
which is much larger than the spin energy. Therefore, 
gradients of x result in rather large increase of total en- 
ergy, and we can neglect them for the trapped system, at 
least for our range of parameters. In other words, healing 
length for \ far exceeds the Thomas-Fermi radius of the 
system, and, therefore, although x is fluctuating inside 
the cloud, it remains nearly constant [24| , except of the 
surface layer, where the density of particles is low. 

Thermal fluctuations of x should be also noticeable in 
three-dimensional condensates at low and moderate tem- 
peratures. In general, the dependences of the number of 
condensed particles on the reduced temperature and crit- 
ical temperature on the total number of atoms for 3D case 
are similar to that in 2D system, which are described by 
Eqs. (2) and (3). The main difference is the powers of t 
and N in the right hand sides of Eqs. (2) and (3) that 
are 3 and —1/3 (hujj_/kT c ~ iV -1 / 3 ), respectively. How- 
ever, in this case one has to take accurately into account 
the possibility of long wave length fluctuations of x i n 
longitudinal direction and formation of kinks |24| . 



IV. CONCLUSIONS 

In this paper, we have studied the effect of thermal 
fluctuations on the coherence among different compo- 
nents of the order parameter in quasi 2D rotating F = 1 
Bose-Einstein condensate, when all three hyperfine states 
are populated. Different axially-symmetric vortex phases 
were considered. We have shown that the deviation of the 
relative phase x = 25o — Si — 5— i from the equilibrium 
can be very significant even at low temperatures, much 
smaller than T c . Fluctuations of relative angle induce siz- 
able fluctuations of the spin texture, namely, local trans- 
verse magnetization of the condensate. We have shown 
that these fluctuations are much more pronounced in an- 
tiferromagnetic case than in the ferromagnetic one. The 
recently proposed in Ref. |20| direct and nondestuctive 
method for the imaging of spinor BEC spatial magnetiza- 
tion (or some of its modification) can be applied for the 
experimental study of the thermal fluctuations of spin 
textures, since it enables multiple-shot imaging and one 
can directly observe the dynamics of a single sample. 
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V. FIGURE CAPTIONS 

Fig. 1. Dependences of Ax (in degrees) on the num- 
ber of particles in the system for different vortex phases 
at fixed value of interaction constant g n — 10 (see in the 
text) and t = 0.1. In the (-1, 0, 1) phase, g s = 0.01g n , 
M = 0.1; in the (1, 1, 1) state, g s = 0.01g n: M = 0.1; in 
the (0, 1, 2) state, g s = — 0.01g n , M = 0.5. Inset shows 
A% as a function of temperature for (0, 1, 2) phase (curve 
1), (1, 1, 1) phase (curve 2), and (-1, 0, 1) phase (curve 
3) at the same values of M, g n and t. Number of atoms 
is N = 1000, interaction constants are g s = — 0.05g n for 
the first curve and g s = 0.05g n for two others. 

Fig. 2. Dependences of (S \ lt r \)x / \hr\ on the total 
magnetization for (-1, 0, 1) phase at N = 1000; g n = 10, 
t = 0.1, g s = 0.01(? n . 



